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3.  Status  of  effort:  A  brief  statement  of  progress  towards  achieving 
the  research  objectives.  (Limit  to  200  words.) 

On  the  theoretical  side,  a  new  highly  scalable  code  for  frozen  core  second  order 
perturbation  theory  gradients  for  closed  shell  molecules  has  been  developed.  The 
development  of  analogous  codes  for  molecules  with  unpaired  electrons  is  in  progress.  A 
working  unrestricted  second  order  perturbation  theory  code  is  now  in  GAMES  S,  and 
work  on  a  scalable  version  of  this  code  will  begin  shortly.  The  derivation  for  the  spin- 
restricted  open  shell  second  order  perturbation  theory  gradients  has  been  completed. 

Other  important  developments  are  new  convergence  methods  for  MCSCF  wavefunctions 
that  facilitate  MCSCF  calculations  on  large  molecules,  the  derivation  of  gradients  for 
multi-reference  second  order  perturbation  theory,  further  developments  of  our  effective 
fragment  potential  (EFP)  method  for  studying  solvation  and  liquid  behavior,  the 
development  of  molecular  dynamics  and  Monte  Carlo  methods  to  facilitate  the  study  of 
solvation  and  liquid  behavior,  the  development  and  implementation  of  a  new  method  for 
producing  global  potential  energy  surfaces  from  sets  of  ab  initio  points,  the  development 
and  implementation  of  a  gridless  approach  to  density  functional  theory,  and  the 
development  and  implementation  of  a  general,  all  electron  MCSCF  approach  to  spin-orbit 
coupling.  With  regard  to  applications,  considerable  progress  has  been  made  in  our 
understanding  of  the  mechanisms  for  formation  of  POSS  (polyhedral  oligomeric 
silsesquioxanes)  reaction.  Careful  studies  of  the  potential  energy  surfaces  of  B  and  B2 
interacting  with  H2  have  been  carried  out,  and  similar  studies  for  A1  have  been  initiated. 
The  reactions  leading  to  A1  and  B  oxides  when  these  species  are  burned  in  a  rocket 
engine  are  underway.  Several  studies  of  high  energy  organic  and  inorganic  compounds 
are  ongoing. 


4.  Accomplishments/New  Findings:  Describe  research  highlights,  their 
significance  to  the  field,  their  relationship  to  the  original  goals,  their  relevance  to 
the  AF's  mission,  and  their  potential  applications  to  AF  and  civilian  technology 
challenges. 

A.  THEORY/MODEL  DEVELOPMENT 

Parallel  Developments.  A  highly  scalable  code  for  frozen  core  second  order  perturbation 
theory  gradients  for  closed  shell  molecules  has  been  developed  and  is  on-line  for  general 
use  at  several  IBM  SP2  and  Cray  T3E  systems,  including  the  T3Es  at  ERDC  and  NAVO 
and  the  SP2  at  Maui.  This  new  code  uses  the  Distributed  Data  Interface  (DDI)  so  that  the 
large  arrays  do  not  have  to  be  replicated.  This  means  that  much  larger  molecules  can  be 
studied  using  geometries  based  on  correlated  wavefunctions.  The  code  scales  very  well 
for  as  many  as  5 12  nodes  and  has  already  been  applied  to  several  challenging 
compounds.  The  development  of  analogous  codes  for  molecules  with  unpaired  electrons, 
using  both  restricted  and  unrestricted  wavefunctions,  is  in  progress.  A  paper  on  the 


restricted  open-shell  derivation  has  been  accepted.  New  developments  will  include  the 
use  of  DDI  for  both  SCF  (completed,  paper  in  print)  and  MCSCF  wavefunctions  (in 
progress).  In  an  important  related  development,  parallel  GAMES S  now  runs  on  clusters 
of  PCs  and  Macs  (running  Linux)  and  high  end  workstations.  This  is  more  challenging 
than  self-contained  massively  parallel  computers,  since  the  overhead  due  to  inter-node 
communications  is  more  complex.  This  issue  has  been  solved  to  some  extent  by  using  a 
Gigabit  Ethernet  switch  with  large  data  packets.  It  is,  however,  still  a  major  research 
area.  As  for  all  other  GAMESS  developments,  we  will  make  our  experiences  in 
developing  scalable  PC  and  Mac  clusters  available  to  all  users. 

Condensed  Phase  Methods.  We  have  already  shown  that  our  effective  fragment  potential 
(EFP)  method  for  solvation  is  excellent  for  water,  in  a  variety  of  applications.  We  are 
now  working  on  extending  the  capabilities  of  the  method  in  several  ways.  We  are 
exploring  several  alternative  approaches  for  incorporating  dispersion  and  other  higher 
order  terms  into  the  method.  Such  terms  are  particularly  important  for  nonpolar  solvents. 
We  are  also  in  the  process  of  extending  the  model  so  that  it  is  equally  applicable  and 
accurate  for  any  solvent.  Key  to  the  success  will  be  the  derivation  of  general  expressions 
for  charge  transfer  and  dispersion  contributions  that  contain  no  fitted  parameters.  A  very 
important  new  development  (by  co-worker  Jan  Jensen)  is  a  new  method  for  using  EFPs 
across  covalent  bonds.  This  will  facilitate  the  representation  of  large  substituents,  as  well 
as  the  treatment  of  large  biological  molecules.  We  have  also  interfaced  the  EFP  method 
with  two  continuum  methods,  the  simple  Onsager  reaction  field  and  the  more 
sophisticated  polarizable  continuum  model  (PCM). 

As  the  number  of  solvent  molecules  in  the  system  increases,  the  number  of 
configurations  to  be  considered  increases  rapidly,  and  traditional  small  molecule 
geometry  optimization  methods  are  not  effective.  We  are  therefore  developing  both 
molecular  dynamics  and  Monte  Carlo  simulation  codes  so  that  the  configurational  space 
can  be  probed  more  effectively,  not  only  for  minima,  but  for  transition  states  and  reaction 
paths  as  well.  This  is  also  an  important  first  step  in  the  development  of  methods  based  on 
the  EFP  model  for  predicting  bulk  properties  and  super-critical  behavior.  We  are 
therefore  incorporating  Ewald  summations  into  our  MD  method.  Collaborators  in  this 
effort  are  Dr.  Paul  Day  (AFRL)  and  Professors  Teresa  Head-Gordon  (UC-Berkeley)  and 
Grant  Merrill  (University  of  Missouri).  Related  to  these  dynamical  methods  is  our 
development  and  implementation  of  a  method  for  converting  the  large  numbers  of  points 
generated  in  ab  initio  trajectories  into  a  global  potential  energy  surface,  using  a  modified 
Shepard  interpolation  approach. 

Major  James  Shoemaker’s  Ph.D  dissertation  in  Engineering  Physics  at  the  Air  Force 
Institute  of  Technology  focused  on  the  development  and  applications  of  an  embedded 
cluster  model  for  treating  surface  chemistry.  A  paper  on  the  theoretical  method,  called 
SIMOMM,  has  been  published  in  the  Journal  of  Physical  Chemistry.  Several  papers 
applying  this  method  have  now  been  published  or  submitted,  and  a  manuscript  that 
describes  applications  to  silicon  carbide  surfaces  is  in  progress.  Extensions  of  this 
method  to  transition  metal  catalysts  are  planned;  however,  molecular  mechanics  is  not 
likely  to  be  a  viable  approach  for  the  bulk,  since  electrons  in  such  systems  are  too 


delocalized.  To  overcome  this  potential  problem,  we  have  begun  to  explore  fast 
multipole  methods  (FMM)  that  scale  linearly  and  are  also  highly  parallel,  so  that  the  bulk 
part  of  the  system  can  be  treated  by  quantum  mechanics.  One  preliminary  paper  has  been 
published  and  the  first  FMM  paper  has  also  been  published. 

Other  developments  There  are  many  problems  that  involve  species  for  which  one  cannot 
write  one  simple  Lewis  structure.  Such  species  must  be  treated  with  multi- 
configurational  (MCSCF)  wavefunctions.  However,  the  computational  demands  for 
MCSCF  calculations  are  very  high.  So,  we  have  developed  a  set  of  methods  for 
converging  such  wavefunctions  that  range  from  very  rapid  convergence  with  high  CPU 
and  disk  requirements  to  methods  that  converge  more  slowly  but  have  very  modest 
computational  demands.  The  advantage  of  the  latter  is  that  much  larger  molecules  than 
normal  can  be  treated  with  this  method. 

For  species  that  must  be  treated  with  MCSCF  wavefunctions,  non-traditional  methods  are 
also  needed  to  recover  the  dynamic  correlation.  One  method  that  has  been  implemented 
into  GAMESS  is  multi-reference  second  order  perturbation  theory.  The  gradient 
formulas  for  this  method  have  been  derived,  and  the  implementation  into  GAMESS  is  in 
progress.  The  key  collaborator  in  this  effort  is  Dr.  Haruyuki  Nakano  (University  of 
Tokyo). 

There  are  many  problems  that  involve  multiple  electronic  states  on  multiple  potential 
energy  surfaces.  Since  such  surfaces  can  cross  and  therefore  induce  radiationless 
transitions,  it  is  important  to  have  the  capability  to  predict  non-adiabatic  coupling  among 
these  surfaces.  One  origin  of  non-adiabatic  interactions  among  surfaces  is  spin-orbit 
coupling.  Spin-orbit  coupling  matrix  elements  must  be  calculated  using  multi-reference 
wavefunctions,  but  until  recently,  there  have  been  no  generally  available  codes  for  the 
computation  of  these  matrix  elements  using  all  electrons  and  any  arbitrary  MCSCF 
wavefunction  for  any  arbitrary  set  of  spin  states.  We  have  now  developed  and 
implemented  such  a  method.  It  will  be  particularly  valuable  for  the  study  of  high  energy 
density  materials  and  for  the  study  of  any  processes  involving  transition  metals,  since 
both  types  of  species  involve  multiple  electronic  states.  Since  such  all-electron 
calculations  are  very  demanding,  two  alternative  approaches  to  this  problem  have  also 
been  develop  and  implemented  into  GAMESS  for  use  with  larger  molecules.  One  of 
these  is  called  the  Zeff  method,  in  which  the  two  electron  term  is  ignored  and  a  fitted 
empirical  parameter  is  introduced  into  the  one-electron  part  of  the  Hamiltonian.  While 
this  is  a  very  efficient  approach,  its  accuracy  is  limited.  However,  there  is  an 
intermediate  method,  in  which  only  those  matrix  elements  connecting  orbitals  in  the 
active  space  with  each  other  are  neglected.  We  have  shown  that  in  this  case  excellent 
agreement  with  the  all  electron  method  is  obtained,  without  any  need  for  fitted 
parameters. 

Density  functional  theory  (DFT)  has  become  extraordinarily  popular  in  recent  years, 
because  it  has  the  same  computational  requirements  as  Hartree-Fock  theory,  but  often  has 
a  level  of  accuracy  that  is  comparable  to  second  order  perturbation  theory.  One  difficulty 
with  DFT  is  that  most  implementations  involve  evaluation  of  integrals  using  numerical 


grids,  and  this  is  inherently  unstable.  We  have  developed  and  implemented  a  gridless 
approach  to  DFT,  by  using  the  resolution  of  the  identity  to  evaluate  the  complicated 
integrals.  This  procedure  requires  the  use  of  auxiliary  basis  sets  to  resolve  the  identity. 
After  extensive  analysis  of  the  convergence  of  several  molecular  properties  as  a  function 
of  basis  set,  we  have  recommended  a  series  of  auxiliary  basis  sets  for  use  with  the  first 
row  in  the  periodic  table.  The  next  step  will  be  to  extend  this  basis  set  development  to 
heavier  elements. 

Another  approach  that  is  very  important  for  large  molecules,  especially  molecules  that 
contain  heavier  elements,  is  the  effective  core  potential  (ECP)  method.  This  method 
replaces  inner  shell  electrons  by  model  potentials,  making  the  calculations  much  more 
efficient.  These  potentials  include  relativistic  corrections  for  the  heavier  elements.  To 
make  these  calculations  more  useful,  we  have  developed  much  more  efficient  methods 
for  obtaining  the  first  and  second  energy  derivatives;  that  is,  the  gradients  and  hessians. 
The  increase  in  speed  is  a  factor  of  5-10,  and  the  method  has  been  implemented  to  take 
advantage  of  parallel  computers. 

As  calculations  on  larger  and  larger  species  are  made  feasible  by  advances  in  parallel 
programs  and  other  innovative  algorithms,  it  becomes  increasingly  important  to  develop 
graphical  methods  fpr  interpreting  the  results  of  complex  calculations.  This  has  led  to  the 
continuing  development  of  MacMolPlt,  a  graphical  interface  for  GAMES  S.  MacMolPlt 
provides  the  ability  to  visualize  (and  rotate  in  real  time)  orbitals,  wavefunctions,  electron 
densities  and  electrostatic  potentials  in  two  and  three  dimensions. 

We  have  also  been  exploring  alternatives  and  extensions  to  the  popular  G2  and  G3 
methods  developed  by  Pople  and  co-workers  for  the  accurate  prediction  of  such 
theormodynamic  properties  as  heats  of  formation,  ionization  potentials  and  electron 
affinities.  One  limitation  of  these  methods  is  that  they  are  applicable  only  to  species  that 
are  adequately  described  by  single  configuration  wavefunctions.  This  eliminates  many 
transition  metal  compounds,  diradicals  and  most  transition  states.  We  have  therefore 
initiated,  in  collaboration  with  the  Radom  group,  the  development  of  multi-reference 
analogs  of  the  G2  and  G3  methods.  The  multi-reference  methods  are  based  on  CASSCF 
wavefunctions,  instead  of  Hartree-Fock,  followed  by  multi-reference  perturbation  theory 
(instead  of  MP2)  and  finally  multi-reference  CL 

B.  APPLICATIONS 

POSS  Compounds.  Polyhedral  oligomeric  silsesquioxanes  (POSS)  have  a  wide  variety  of 
important  applications,  including  lubricants  and  coatings.  Consequently,  there  is 
considerable  interest  in  these  species  in  industry,  universities  and  Air  Force  Laboratories. 
We  have  therefore  initiated  a  detailed  study  of  POSS  compounds,  with  primary  emphasis 
on  the  mechanisms  by  which  these  species  form,  in  both  the  gas  phase,  and  in  solution 
and  the  presence  of  catalysts.  In  the  first  series  of  calculations,  we  examined  the  three 
stepwise  hydrolysis  steps  of  trichlorosilane,  followed  by  the  condensation  of  the  product 
trihydroxysilane.  All  four  reactions  are  predicted  to  have  large  barrier  heights  in  the  gas 
phase,  but  the  addition  of  just  one  extra  water  molecule  is  sufficient  to  reduce  all  barriers 


to  zero  or  nearly  zero,  except  for  the  first  hydrolysis  step.  In  the  latter,  there  is  still  a 
residual  barrier  of  a  few  kcal/mol.  So,  the  next  step  will  be  to  examine  the  effect  of 
additional  water  molecules,  especially  on  this  first  step.  Since  these  are  rather  large 
species,  the  additional  water  molecules  will  be  treated  with  the  EFP  method.  The 
molecular  dynamics  (MD)  and  Monte  Carlo  (MC)  codes,  now  almost  completed,  will  be 
essential  as  the  number  of  water  molecules  is  increased.  A  new  paper  discussing  the 
effects  of  waters  and  subituent  effects  is  in  preparation. 

For  both  the  hydrolysis  and  condensation  steps  in  the  POSS  mechanism,  there  exist  low 
lying  intermediates.  Even  though  the  net  barriers  are  nearly  all  small  in  the  presence  of 
an  additional  water  molecule,  the  barriers  separating  these  intermediates  from  products 
are  large.  So,  it  is  not  entirely  clear  how  easily  products  will  be  obtained  in  solvent.  To 
determine  this,  ab  initio  dynamics  “on-the-fly”  will  be  carried  out  in  the  presence  of 
waters,  using  the  MD  and  MC  codes. 

Calculations  have  also  been  completed  on  the  next  steps  in  the  mechanism,  in  which  the 
initial  condensation  products,  disiloxanes,  are  further  condensed  to  the  ring  compounds 
D3  and  D4.  These  results  are  consistent  with  those  reported  above,  in  that  initially  high 
barriers  are  reduced  to  nearly  zero  by  an  additional  water  molecule.  The  next  steps  in 
this  process  will  include  further  condensation  to  the  three-dimensional  cages  and  analysis 
of  substituent  effects  on  the  reaction  mechanism  and  energetics. 

Although  our  main  focus  is  on  Si-based  POSS,  the  Ti  analogs  are  also  of  interest.  We 
have  therefore  initiated  a  series  of  calculations  on  the  Ti-POSS  compounds  to  study  their 
properties. 

Fuels  with  High  N-Content.  There  has  been  considerable  interest  in  the  last  several  years 
in  the  potential  for  cubic  molecules  as  high  energy  fuels.  Cubane  itself  has  been 
considered,  and  so  far  rejected  because  of  the  complexity  of  its  synthesis.  The  octasila 
analog  of  cubane  has  been  prepared.  Unlike  the  carbon  system,  the  cubic  structure  is  the 
most  stable  SigHg  isomer.  This  makes  it  potentially  very  interesting,  since  it  is  still  rather 
high  in  energy.  However,  the  heavy  mass  of  Si  precludes  any  viability  of  such  species  as 
high  energy  fuels. 

On  the  other  hand,  cubic  Ng  is  a  light,  very  high  energy  fuel  that  has  been  the  subject  of 
considerable  attention.  Extensive  calculations,  however,  illustrate  that  the  upper  limit  of 
the  barrier  separating  this  species  from  four  nitrogen  molecules  is  20  kcal/mol.  It  is, 
therefore,  not  a  viable  high  energy  fuel.  In  the  process  of  studying  the  cubic  Ng  potential 
energy  surface,  it  was  discovered  that  there  are  three  Ng  isomers  that  are  much  lower  in 
energy  than  the  cube,  yet  still  much  higher  in  energy  than  molecular  nitrogen.  The 
potential  energy  surfaces  for  these  isomers  have  now  been  studied  in  detail.  Two  of  these 
have  barriers  on  the  order  of  15-20  kcal/mol  (too  small),  but  the  third  has  a  barrier  of 
nearly  30  kcal/mol.  So,  this  isomer  may  be  a  viable  synthetic  target. 

Related  to  the  foregoing,  much  attention  has  been  paid  recently  to  N4  as  a  potential  fuel. 
Many  references  are  made  to  a  barrier  on  the  order  of  50-65  kcal/mol,  even  though  it  has 


been  shown  by  Yarkony  that  non-adiabatic  interactions  reduce  this  barrier  to  less  than  40 
kcal/mol.  For  the  first  time,  we  have  mapped  out  the  potential  energy  surface  that  leads 
from  tetrahedral  N4,  using  second  order  perturbation  theory.  One  also  must  consider  the 
possibility  that  two  N4  molecules  could  interact  and  form  molecular  nitrogen.  This 
possibility  is  being  explored  in  this  laboratory. 

Rob  Schmitt  at  SRI  proposed  an  intriguing  series  of  potential  high  energy  compounds 
that  we  have  been  exploring  computationally.  The  procedure  is  to  use  a  modified  version 
of  G2  to  predict  the  heats  of  formation  and  then  to  use  these  predictions  to  assess  the 
specific  impulse.  For  the  most  promising  species,  further  studies  are  performed  in  which 
the  stability  of  these  compounds  with  respect  to  unimolecular  decomposition, 
bimolecular  interactions  and  attack  by  environmental  species  such  as  water. 

Christe  and  co-workers  at  AFRL  (Edwards)  have  isolated  the  first  new  all-nitrogen 
species  in  nearly  a  century.  With  the  counter  ion,  this  species  is  [AsF6'][N5+].  We  have 
now  studied  the  potential  energy  surface  of  this  species  as  it  forms  AsFs  +  FN5,  and 

subsequently  the  pathway  leading  from  FN5  to  FN3  +  N2. 

There  has  been  considerable  experimental  and  theoretical  interest  in  metal-doped  solid 
hydrogen.  We  have,  an  ongoing  series  of  investigations  of  the  low-lying  electronic  states 
of  BH2  system,  to  determine  the  energetics  by  which  the  metal  is  held  as  a  weakly  bound 
species  in  the  matrix  and  the  ease  with  which  two  metal  atoms  find  each  other.  Similar 
studies  are  under  way  on  the  B2-H2  system  and  the  A1  analogs. 

The  desired  reaction  product  for  burning  Al-doped  H2  is  AI2O3,  since  it  is  the 
thermodynamically  most  stable  aluminum  oxide.  However,  it  is  unknown  how  or  when 
this  species  forms  in  the  process.  We  have  therefore  initiated  an  extensive  study  of  the 
potential  energy  surface  that  leads  from  A1  and  02  to  various  oxides.  The  first  step,  the 
reaction  of  A1  with  02  to  from  A102  is  nearly  complete.  The  disproportionation  reaction 
between  two  A102  molecules  is  in  progress. 
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